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Abstract 

QED processes at electron-positron colliders are considerd. We present differential 
cross-sections for large-angle Bhabha scattering, annihilation into muons and photons. 
Radiative corrections in the first order are taken into account exactly. Leading loga- 
rithmic contributions are calculated in all orders by means of the structure-function 
method. An accuracy of the calculation can be estimated about 0.2%. 

PACS 12.20.-m Quantum electrodynamics, 12.20.Ds Specific calculations 

1 Introduction 

At the existing VEPP-2M e + e~ collider and planned meson factories QED processes of the 
lowest order of the perturbation theory (PT) play an important role. These processes are 
to be considered as an essential background when extracting subtle mesons properties from 
experimental data. QED processes are used also for calibration and monitoring purposes. 
For instance, large-angle Bhabha scattering is used for a precise determination of luminosity 
at e + e~ colliders. That is the reason, why radiative corrections (RC) to QED processes are 
to be considered in detail. 

A lot of attention was paid to the problem of radiative correction calculations to various 
QED processes, starting more than 50 years ago ffl, concerning mainly the lowest order of 



PT calculations. The accuracy requirements of modern experiments, however, exceed the 
ones provided by the first order RC. Unfortunately, calculations of radiative corrections in 
higher orders encounter tremendous technical difficulties. Nevertheless, powerful methods 
developed in quantum chromodynamics (see paper [f2| and references therein) provide a 
possibility to improve essentially the results obtained earlier. At first we mean the methods 
based on the renormalization group ideas and on the factorization theorem. Thay allow us 
to get a differential cross-section of a certain process similar to the Drell-Yan process cross- 
section, and to consider leading logarithmic RC of higher orders. Meanwhile, nonleading 
contributions are to be taken from the lowest order PT calculations. 

The aim of our paper is to provide relevant theoretical formulae for QED processes, which 
are required for experiments with CMD-2 and SND detectors at VEPP-2M (Novosibirsk) 
collider, at DA$NE (Frascati) Q and BEPC/BES (Beijing) || machines. Formulae cited 
below may be applied also for higher energies (see the Conclusions), if one will take into 
account weak interaction and higher hadronic resonance contributions. 

A cross-section calculated with an account of radiative corrections (RC) in the n-th order 
of perturbation theory (PT) contains enhanced contributions of the form (a/7r) n L n , where 
L = ln(s/?7ig) is the large logarithm (for s ~ 1 GeV 2 , L « 15). We call these contributions 
the leading ones. Nonleading contributions have the order (a/7r) n L m , m < n. The terms, 
proportional to (a/7r) n L n , can be calculated by means of the structure function method ||. 
The structure function formalism, based on the renormalization group approach, permits 
one to keep all leading terms of order (aL/ii) n , n = 0, 1, 2 . . . explicitly. 

In the first order of PT the nonleading terms, proportional to (a /it), can be accounted 
by means of so-called /C-factor, /C = 1 + (a/ir)K. As for nonleading terms of the order 
(a/ir) 2 L, they can be correctly calculated in the two-loop approximation. We will consider 
them in further publications 0. Fortunately, this second order nonleading RC are small: 
(a/ir) 2 L ~ 10~ 4 . So, the presented formulae guarantee theoretical precision of the order 
0.2%. Considering hard photon emission in the first order of PT we distinguish the cont- 
ributions due to the initial state radiation, the final state radiation, and their interference. 
The first one always contains large logarithms. The ones due to the final state radiation as 
well as the initial-final state interference do not enhance the corrections considerably, except 
the case of Bhabha scattering. In the case of hard photon radiation large logarithms come 
from kinematical regions, where the photona are emitted along electron and positron beams. 
We will call this kinematics as the collinear one. 

Our work consist in explicit calculations of the Born and the first order RC contributi- 
ons to differential cross-sections. We apply the structure function method to increase the 
accuracy. 

The paper is organized as follows. In Section 2 we consider /i + \x~ production. Both 
charge-even and charge-odd contributions to the differential cross-section are evaluated. In 
Sect. 3 we consider the electron-positron scattering process at large angles. In Sect. 4 we 
investigate electron-positron annihilation into photons. In the Conclusions we discuss the 
results obtained and estimate the provided precision. Numerical illustrations are given for a 
realistic experimental set up. 



2 Muon Pair Production 

Consider the process 

e + 0+) + e~0_) -► /i + (g+) + A/"(g_). 

Taking into account the photon and Z-boson intermediate states, the differential cross- 
section in the Born approximation (in the framework of the Standard Model) has the form: 
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where the centre-of-mass reference system of the initial beams is applied. 

Here Ky/ (we put the explicit expression for it in the Conclusions) represents contribu- 
tions due to Z-boson intermediate states (see || for example), we will neglect them within 
the accepted precision: Kw ~ s /^z ~ 10~ 3 . We will drop also all other contributions due 
to weak interactions in higher orders. 

Consider first the even, with respect to the c ^ —c permutation, part of the one-loop 
virtual and soft radiative corrections. Using the known in the literature Dirac and Pauli 
muon form factors and the known soft photon contributions, we obtain 
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where Ae <C £ is the maximum energy of soft photon in the centre-of-mass system. Il(s) is 
the vacuum polarization operator including electron, muon, tau-meson and hadron contri- 
butions (see Appendix 1). 

The odd part of the one-loop correction comes from the interference of Born and box 
Feynman diagrams and from the interference part of the soft photon emission contribution. 
It causes the charge asymmetry of the process: 



rj = [do-(c) - da(-c)]/[da(c) + dcr(-c)] ^ 0. 
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The odd part of the differential cross-section has the following form: 
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In the ultra-relativistic limit {(3 — > 1) for the even part we obtain: 
f <fr£f +v \ = dap 1 
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For the odd part in this limit we obtained the same as Khriplovich [I 
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Consider now the process of hard photon emission 

e + (p + ) + e~(p_) ->■ /J + (g+) + ju~(g_) + j(k). 
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It was studied in detail ||13|| . The photon energy is assumed to be larger than Ae. The 
differential cross-section has the form 
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The contribution due to the initial state radiation reads 
R e = A e + r> e , 
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The final state radiation and the interference of the initial and final state radiation contri- 
butions are 
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Quantity i? e contains collinear and infrared singularities. A^ and A e/1 have only infrared 
singularities. B^ and I? eM are free from singularities. Quantity R in the ultra-relativistic 
case is given in Appendix 3. 

After algebraic transformations we get 
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Note that from these expressions one may in a moment obtain the corresponding matrix 
element of the cross symmetrical process e~u + — > e _ /i + 7. To be rigorous, we have note 
that in the cross symmetrical channel one has to take into account some additional terms 
proportional to m 2 . In our channel they can be shown as follows: 
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We checked the matrix element by a comparison with the one used in FORTRAN program [13 1 . 
describing electron-muon scattering. 

The sum of the hard photon contribution, integrated over the photon phase volume with 
the condition k° > As, and the contribution due to the soft and virtual photon emission does 
not depend on the auxiliary parameter A = Ae/e <C 1. The main contribution, proportional 
to the large logarithm, comes from the integration of R e in the case of collinear kinematics 
of photon emission. For defmiteness let us consider the case when the photon moves close 
to the initial electron direction: 
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where x is the energy fraction carried away by the emitted photon, x = k°/e = 1 — Si/s. 
Performing the integration over the photon emission angles, we can present the corresponding 
part of the cross-section (a similar contribution of the hard photon emission along the 
positron is included below also) in the form 
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where d<7o(l — xi,l — x 2 )/dQ_ is the so-called shifted Born differential cross-section. It 
describes the process e + {p + {l — x 2 )) + e~(p_(l — Xi)) — ► ji + (q + ) + n~(q), 
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Using the conservation laws 

*i + z 2 = Fi + Y 2 , 
y x \j\-c 2 _ = y 2 s]\ - d . 

we obtain the energy fraction of the created muon 
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Quantity C, after adding the corrections due to soft and virtual photons, turns out to be 
the lowest order perturbative expansion of the convolution of the structure function T> with 
the shifted Born differential cross-section. Quantity D plays role of a compensating term. 
Namely, in the sum with the contribution of the cross-section due to hard (A; > As) photon 
emission at angles (with respect to the electron and positron) larger than 9q the dependence 
on the auxiliary parameters will cancel. 

Here we remind about experimental conditions of the final particles detection mentioned 
above. They are to be imposed explicitly by introducing the restriction of the following kind: 
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where y t h£ = £th is the threshold of the detectors, angle \l/o determines the dead cones, 
surrounding beam axes, unattainable for detection. More detailed cuts can be implemented 
in a Monte Carlo program, using the formulae given above. 

The leading contributions to the cross-section, containing large logarithm L, as may be 
recognized, combine to the kernel of Altarelli-Parisi-Lipatov evolution equation: 
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This formula is valid in the leading logarithmical approximation. We will modify it by in- 
cluding nonleading contributions and using the smoothed representation for structure func- 
tions HI: 
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In comparison with the corresponding formula in ref. we shifted the terms, arising due to 
virtual e + e~ pair production corrections, from D 7 into T> e e . 

Finally, the differential cross-section can be presented in the form 
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where D, C e/1 and C M are compensating terms, which provide cancellation of auxiliary pa- 
rameters A and 9q inside figure brackets. In the first term, containing T> functions, we 
gather all leading terms. A part of nonleading terms proportional to the Born cross-section 
is written as the /C-factor. The rest nonleading terms are written as two additional terms. 
The compensating term D (see Eq. ([16])) comes from the integration in the collinear region 
of hard photon emission. Quantities C^ and C e ^ come from the even and odd parts of the 
differential cross-section (arising due to soft and virtual corrections), respectively. Here we 
consider the phase volumes of two (df2_) and three (dr) final particles as the ones, which 
already include all required experimental cuts. Using the conservation laws Eq. fll8|) and 
concrete experimental conditions one can define the lower limits of the integration over z\ 
and z 2 . 

There is a peculiar feature in the spectrum of hard photons. Namely, in the end of the 
spectrum the differential cross-section is proportional to the factor 
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which defines a peak at Sj « 5.6m 2 It comes from the Feynman diagrams describing the 
emission by the initial particles fll5"|. 



3 Large— Angle Bhabha Scattering 

The cross-section of Bhabha scattering (corrected by the vacuum polarization factor), which 
enter into the Drell-Yan form of corrected cross-section, has a bit more complicated form, as 
far as the scattering and annihilation amplitudes and their interference are to be taken into 
account. We remind here the form of the Lorentz-invariant matrix element module squared 
in the Born approximation: 
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The first term in the right hand side describes the scattering-type Feynman diagram square. 
The second one corresponds to the square of the annihilation-type diagram. And the third 
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one deals with the interference of the two diagrams. A more compact representation of Rq is 



also useful, Ro 
has the form 



s/t + t/s) 2 . The differential cross-section in the Born approximation 
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We will need also quantity R for arbitrary energies of initial particles. Suppose that the 
initial electron and positron lost a certain energy fraction. The corresponding kinematics is 
defined as follows: 
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Here the shifted Born cross-section corrected by vacuum polarization insertions into virtual 
photon propagators reads 
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Rewriting the known results ]TT] , [T2| for the cross-section in the Born approximation 
with one-loop virtual corrections to it and with the other ones arising due to soft photon 
emission, we obtain 
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e, 
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is the part of the /C-factor coming from soft and virtual photon corrections, 
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Quantity Ae in Eq. (|2~8|) is the maximum energy of emitted soft photons. II(s) and II(t) 
are the vacuum polarization operators in the s and t channels. In the Conclusions we will 
estimate the contribution of weak interactions. 

Consider now the process of hard photon (with the energy u = k° > Ae) emission 

e + (p+) + e~(p_) -> e + (p' + ) + e~(p'_) + ^(k). 



We start with the differential cross-section in the form suggested by F.A. Berends et al. |TT 
(which is valid for scattering angles being large compared with m e /e): 

1 ^ P AT AT d 3 p' + d 3 P '_d 3 k (4) , , 

dcxhard = tt^ - R en dT, dT = , , — b K >{p+ + p--p + -p_- k), (31) 

WT m\ (s t \ 2 ml ( s ti x2 

^ " ^"(^U + ; + 1 J ~(^U + 7 + 1 

m 2 ( S\ t \ 2 m 2 { si t 1 ' 



where 






W = 1 1 1 — 

X+X- X'+X- X'+X+ X'-X- X'+X- X-X+ 



T 



ssi(s 2 + s\) + tti(t 2 + t\) + uuxiu 2 + u\) 



SSitti 



and the invariants are defined as 



s = 2p_p+, si = 2pV+, t = -2p_pi, t\ = -2p + p' + , 
u = -2p_p' + , m = -2p + p'_, x± = kp±, X± = kp'±. 

It is convenient to extract the contribution of the collinear kinematics. We do that for 
the following reasons. First, it is natural to separate the region with very a sharp behaviour 
of the cross-section and to consider it carefully. Second, we keep in mind the idea of the 
leading logarithm factorization, which is valid in all orders of the perturbation theory. We 
will evaluate the collinear kinematical regions in two different ways. The first one (the 
quasireal electron approximation) is suitable for a generalization in order to account higher 
order leading corrections by means of the structure function method. In this way we will 
obtain below the leading logarithmic contributions and the compensating terms, which will 
provide the cancellation of auxiliary parameters. The second one (the direct calculation) is 
more rigorous, it can be used as a check of the first one. We discuss it in detail in Appendix 2. 

To obtain explicit formulae for compensators it is needed to consider four kinematical 
regions corresponding to hard photon emission inside narrow cones, surrounding the initial 
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and final charged particle momenta. The vertices of the cones are taken in the interaction 
point. We introduce a small auxiliary parameter 9 , it should obey the restriction 



m, 



s <6> < 1. 



(32) 



So, we define a collinear kinematical region, as the part of the whole phase space, in which 
the hard photon is emitted within the cone of 8q polar angle with respect to the direction of 
motion of one of the charged particles. 



Using the method of quasireal electrons |16j , the matrix element A4 (squared and summed 
up over polarization states) of the process of hard photon emission can be expressed through 
a shifted matrix element of the process without photon emission (see Eq. (4) in |TjJ ) : 



J2\M( Pu k,p' v X)\ 2 
52\M(pi,A,k,X)\ a 



4ira 
4ira 



x\ 



m 



x(l — x) kpi (kpi) 2 



£|M,(Pi-*,Pi,*)| ! 



y 2 + Y 2 e 



m 



x 



loY kpi (kpi) 2 
o 
:. a —, Y = x + y 



J2\M (p 1 ,p[ + k,X)\ 2 , (33) 



w o Pi 

-, Pi = e, y = — 

e e 

where X denotes the momenta of non-radiating incoming and outgoing particles in a concrete 
process. The integration over the phase volume of the emitted photon inside the narrow cone, 
surrounding its parent charged particle momentum, gives the following factors: 



4a f d 3 k 

167T 2 J CO 
+ 1-Z1 

4a f d 3 k 

167T 2 J CO 
-1-9 In 7o 


"1 + (1-: 

x(l — X 

, Zi = 1 

y + y 2 


c) 2 1 m 2 


a dz\ 


) fcpi (h 

-x, 

1 m 2 


l )\ 


2tt 
a d^3 


Z\ 

p(l 

= 1 


xY I 


zp'i (kp[) 2 

Z3 = l 


2tt z 3 

CO 


-|-Z 111 23 ti ^3 


p\ 


° + co 



Pg\ Zl )(L-l-\ In- 



:, 3 )(^-i+in| 

x 



9, 



(34) 



Y' 

Note that the terms proportional to (L — 1) contain the kernel P^ (see Eq. fl21"D) of Altarelli- 
Parisi-Lipatov evolution equations (more precisely, they contain 0-part of the nonsinglet 
kernel): 

1 + z 2 

-G(l-z-A). 



Pi l \z) 



1- z 



Collecting the contributions of the four collinear regions, we obtain 



dcr, 



coll 



dfi_ 



a r dx 



IT j x 
A 



1 " x + t) ( l " 1 + ln f + 2 ln(1 " x) ) + y 



d*o(l,l) 

dcT (l — x, 1) dcr (l, 1 — x) 



\-,+*L-^i) + * 



dQ. 



dfi_ 



(35) 
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where the shifted Born cross-section is defined in Eq. Q2"7|). 

Adding the contributions of virtual and soft photon emission, we restore the complete 
kernel. Generalizing the procedure for the case of photon emission by all charged particles, 
we come to the representation of the cross-section in the leading logarithmic approximation. 
The final expression for the cross-section therefore has the form 



da e+e ~* e+e W 



Yi 



1 1 

J 'dzi fdz 2 V(z 1 )V(z 2 



da (z h z 2 ) f 1+ a Ksv]() 



dfi_ 



n 



2/th 2/th 




ar . , 
x H I In 



^ 2 fl - aO 



ln^ + 

4 



.r 



2 n 



a? 



4« 2 



dcr, 



Born 




dfi_ 



s(l-x) 2 [2-x(l-c)] 4 



'3 - 3a; + x 2 + 2a;(2 - x)c + c 2 (l - x + x 2 



4a* 



li 



2-1 



' s[2- 


-x(l + c 


■)}'\ 


a 2 
4s 


2z x z 2 


/ 7T 


2l + 2 2 - c(zi 

2/th(l + c) 


-2 2 ) 



1-c / 

3 - 3a; + x 2 - 2x(2 - x)c + c 2 (l - x + x 2 )\ 2 

ln(ctg-)ln— + ^-^- / — r- 0- 



>0 



^ 



2vr 2 s y 4 

fc°>Ae 
7r-e ( ,>6»>6» 

2? + 4 - ( z i - 4) C 



da 



(36) 



2 - y th (l - c) 



2 2 



2l + 2 2 - C(Zl - Z 2 J 

Zi!/th(l - c) 



22i - y t h(l + c) 



The last term describes hard photon emission process, provided that the photon energy 
fraction x is larger than A = Ae/e, and its emission angle with respect to any charged 
particle direction is larger than some small quantity 6 . The sum of the last 3 terms in 
Eq. ( |36|) does not depend on the auxiliary parameters A and 9q, if they are sufficiently small. 
We omitted the effects due to vacuum polarization in the last three terms which describe real 
hard photon emission. Because the theoretical uncertainty, coming from this approximation, 
has the order 5{da)/da ~ (a/n) 2 L < 1CT 4 . Nevertheless if the centre-off-mass energy is 
close to some resonance mass (say to m^) the effect due to vacuum polarization may become 
visible. The differential cross-section of non-collinear hard photon emission, that takes into 
account vacuum polarization explicitly, is presented in Appendix 3. 
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4 Annihilation of e + e into photons 

Considering the RC due to emission of virtual and soft real photons to the cross-section of 
two quantum annihilation process 



e + (p + ) + e (p_ 
we will use the results obtained in papers Jl 

a 



7(9i) + 7fe), 



d&B+S+V 

K S v 

da (l,l) 



d<7 (l,l) 1 + 



7T 



ZL 1 ~ c 

y + 2(i+c 2 ) 






(37) 



(38) 



2 1 



+ 1 + 



1 



c 



1 + c 2 1 
a 2 (l+c 2 ) 



11+ c\ 2 1 - c 

In + (c 



s(l 



dfii , s = (p + +p_) , c = cos6'i, 61 = (fj>_ 



We suppose that the two final photons are registered in an experiment and their polar angles 
with respect to the initial beam directions are not small (#1.2 3> m e /e). 
Consider the three-quantum annihilation process 

e + (p + ) + e~(p_) -► 7(gi) + 7(^2) + 7(93) 
with the cross-section (see the paper by M.V. Terentjev |L] 



da' 



>3 7 



«" 



R 



3o 



dr 



8 X§ + (x 3 ) 2 

X1X2X1X2 



(39) 



2mi 



xl + xl 



/ \2 



-XiX2(Xs) 2 
two cyclic permutations, 

d 3 gid 3 g 2 d 3 g 3 



+ 



(x'i) 2 + (x' 2 ) 
X1X2X3 



9l92?3 



-<5 (4) (p++p_ -gi -g 2 -g 3 ), 



where 



Xi = 9iP- 



Xi = ?#+, 



1,2,3. 



The process can be treated as a radiative correction to the two-quantum annihilation. 

In the same way as we have done before we will distinguish the contributions of the 
collinear kinematical region, when extra photons are emitted within narrow cones of the 
opening angle 29$ < 1 to one of the charged particles and the semi-collinear ones, when 
extra photons are emitted outside these cones. This contribution can be obtained using the 
quasireal electron method [Ol. It reads: 



do-coii 



a 

n 



dx 
x 



[1-X + 



X 



L 



4 



+ 



X 



(40) 



x [da (l-x, l)+dcr (l, 1-x)], 
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where the shifted cross-section has the form 

la 2 z\(\ - cf + 4(1 + c ) 2 



da (z 1 ,z 2 ) 



s (l-c 2 )(z 1 + z 2 + (z 2 -z x )c) 2 



dfli 



(41) 



Again rearranging the separate contributions and applying the structure functions method, 
we obtain the improved cross-section 



da e 



+77 ( 7 ) 



a r dx 



dzi V(zi) J dz 2 V(z 2 )da (z 1} z 2 ) ( 1 + -K sv 



a. 



21 



22 

„2\ /)2 



Z; 



+ 

IT J X 
A 

1 

+ 3 

z t >A 

TT-8 o >e t >0 o 



x"\. 0Z x' 

l-x-\ In — H 

2 ) 4 2 



do"o(l — x, 1) + d<7o(l, 1 — x) 



4a 3 

TT 2 S 2 



4(i + 4) 



zlzl{l-cl){l-cl 



two cyclic permutations 



dr, 



(42) 



cos6»i, 9i=p_q. 



where lower limits z± j2 are defined in Eq. (|36D . The multiplier | in the last term takes into 
account the identity of the final photons. The sum of the last two terms does not depend 
on A and #o- Note that the annihilation process is a pure QED one, hadronic contributions 
as well as weak interaction effects are far beyond the required accuracy. 



5 Conclusions 

Thus we had considered the series of processes at electron-positron colliders of moderately 
high energies. We presented differential cross-sections to be integrated over concrete expe- 
rimental conditions. The formulae are good as for semi-analytical integration, as well as 
for the creation of a Monte Carlo event generator ||20|| . In a separate publication we are 
going to present analysis of the effects of radiative corrections for the conditions of VEPP- 
2M (Novosibirsk), DA$NE (Frascati) and BEPC/BES (Beijing). The idea of our approach 
was to separate the contributions due to 2 — > 2 like processes and 2 — ^ 3 like ones. The 
compensating terms allow us to eliminate the dependence on auxiliary parameters in both 
contributions separately. 

Note that all presented formulae are valid only for large-angle processes. Indeed, in the 
region of very small angles 6 ~ m e /e of final particles with respect to the beam directions 
there are contributions of double logarithmic approximation [9] . These small angle regions 
give the main part of the total cross-section. We suppose that this kinematics is rejected by 
experimental cuts. 

In the Table and Figures we present some results of numerical calculations according to 
our formulae. We suppose that a process-event implies detecting of two final particles with 
the polar angles with respect to the beam axes more than some value \l/o- The energies of 
the particles have to exeed some experimental threshold e t h- A cut-off on the acollinearity 
of the final particle momenta is possible. But we switched off it in the computations. 
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Table 1: The values of Bhabha cross-section and radiative corrections to it. 



e± 


a$ OTn (mb) 


Svp (%) 


S& (%) 


*SF (%) 


5k (%) 


5 1 (%) 


cr to t (mb) 


9° < 6± < 171° 


3.77- 1(T 2 


1.53 


0.11 


-0.30 


-1.82 


1.50 


3.81 ■ 10~ 2 


1° < 6± < 179° 


3.08 


0.80 


0.70 


-0.17 


-3.5 


4.74 


3.16 



In the Table we give the values of different RC contributions to large-angle Bhabha 
scattering cross-section. Here we switched off also the cut-off on the final particle energies 
(we used only kinematical restrictions). The contributions (see Eq.fl36|)) are defined as 
follows: 



Vtot 



da e 



"(7) 



dfi_ 



da 



J Born 

^°— da 
da 



i + 



100% 



VP 



^SF + ^SF + ^K + & 



,(43) 



where £yp is due to the vacuum polarization being included into the Born level diagrams; 
<5gp is due to the initial (final) state leading logarithmic corrections; 8k shows the impact 
of the /C-f actor; 8 1 describes the contribution of one hard photon emission at large angles. 
The effect due to the width of meson is included as a part of vacuum polarization. It 
is small for the given integrated cross-sections. But for the description of a differential 
cross-section it is important, especially for large scattering angles (see |2c 



M) 



In Figures 1 and 2 we illustrate the charge-odd part of the differential cross-section for 



the e + e — > fi + fi process (see Eq.(|3|)). The quantity 



A 



da. 



+/X+/U (7) 



FB 



odd 



/dc 



da, 



e+e — —rfj,+/i~ (7) 




/dc 



100% 



(44) 



is shown there as a function of c. The short-dashed line represents the contribution due to 
soft photon emission and virtual corrections. The long-dashed line represents the corres- 
ponding odd contribution due to hard photon emission. It comes from the intereference of 
the amplitudes due to initial and final radiation (R efl ). In the sum of the two contributions 
the dependence on the auxiliary parameter A disappears (it was chosen A = 0.01). And we 
obtain an experimentally measurable asymmetry Aps (the solid line). There is also a con- 
tribution to the asymmetry due to electroweak interactions. Namely, due to the interference 
of the Born level amplitudes with 7 and Z boson in the s-channel. It can be found from 
the weak /C-factor fl46|) . It is included in the total sums (solid lines). But for the chosen 
energies it is really small (it gives a maximal shift of about 0.01% for i?beam = 0.51 GeV and 
about 0.1% for £ beam = 1.55 GeV. 

Let us discuss the accuracy, provided by our formulae. The contribution of weak interac- 
tion to the cross-section of muon pair production and Bhabha scattering was parameterized 
be so-called weak /C-factors: 



K. 



w 



(d<r) E w - (do-) QED 



(da) 



QED 



(45) 



where quantities (dcr)Ew and (da)QED are the cross-sections calculated in the Born appro- 
ximation in the frames of the Standard Model and QED, respectively. The weak /C-factors 
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(we used the results of papers 



arc 






K W 



s 2 (2-/3 2 (l-c 2 ))- 1 f, „,, .././ M 2 Z \ 

)(2_/3 2 (l-c 2 )) lc 2 v (3-2^)+c 2 a 



(s - Ml) 2 + Mfr| 
^(cl + cD + cP 



lll-Mu+8^ 



2 sin 2$^ 
(1 - c) 2 



2(3 + c 



2^|2 



c v = c (l -4sin 2 6> w ), 



45j + (1 - c) 2 5 2 + (1 + c) 2 £ 3 



-1, 



(46) 



(47) 



B 3 
X 
A 



i + {gl - gl)i 



Bo 



As 



1 + (9v ~ 9 2 a)X 



1 + S ~ + (9v - 9af(- t i + x) 



M 



s — m 2 + iM z T z 
G F M 2 Z 



t-M 2 ' 



2^/2 



(sin 29 



w) 



na 



9a 



2 



9v 



(l-4sin 2 w ), 



here 9y/ is the weak mixing angle. 

These quantities are of order 0.1% up to *Js < 3 GeV. Contribution of weak interactions 
to the cross-section of the annihilation into photons (which is absent at the Born level) can 
be estimated as 



(K 



W )ee— >77 



< 



as 



(4* 



It comes from one-loop electroweak radiative corrections. Another source of uncertainties 
comes from the approximation of collinear kinematics (or the approximation of quasireal 



electrons f[Uj) It can be estimated by the largest omitted terms 



*9l and «(!£ 

7T \e9r) 



7T 



(49) 



Really in the calculations we used the value of the parameter #o of the order 10~ 2 because 
of the restrictions 8 ^ 1 and s6o/m e ^> 1. Note that the coefficient before terms of that 
sort (calculable in principle) is the function of energy fractions and angles, they are of order 
of 1. For typical values of energy e = 0.5 GeV this uncertainty is of order 2 ■ 10~ 4 or less. 

The third source of uncertainties is the error in the definition of the hadronic vacuum 
polarization. It has been estimated || to be of order 0.04%. For 0-meson factories a 
systematic error in the definition of the 0-meson contribution into vacuum polarization is 
to be added. 

Next point concerns nonleading terms of order (a/ir) 2 L. There are several sources of 
them. One is the emission of two extra hard particles (for the case of Bhabha scattering it 
was considered in the series of papers |ll[). Other are related to virtual and soft-photon 
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radiative corrections to single hard photon emission and Born processes. The most part 
of these contributions was not considered up to now. Nevertheless, we can estimate the 
coefficient before the quantity (a/ir) 2 L p» 1 • 10~ 4 to be of order of unity. That was indirectly 
confirmed by our complete calculations of these terms for the case of small-angle Bhabha 
scattering [^1 . 



Considering all mentioned above sources of uncertainties as independent, we conclude 
that the systematic error of our formulae does not exceed 0.2%. The main error is due to 
unknown second-order next-to-leading radiative corrections. 

For precise luminosity measurements we suggest to use the large-angle Bhabha scattering 
process. It has a very large cross-section, a good signature in detectors, and the lowest 
theoretical uncertainty. 
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Appendix 1 

We present here leptonic and hadronic contributions into the vacuum polarization operator: 



n 
n*00 

U h (s) 



-n 1 ( s ) + f-) 2 n 2 ( s ) + f-) 3 n 3 ( s ) + ... 

7T \7T/ \TT / 



(50) 



Ana 



PV 



^.e+e — >hadrons / „/"\ 



-ds' — ma 



e + e — +hadrons, 



iml 



S — S 



s) 



The first order leptonic contribution is well known |T|: 

1 



ni( s 

where 



1 5 

^L - - + f{x„) + f{x T ) - ivr 



+ 0(xj9(l - X,) + (J)(X T )Q(1 - X T ) 



(51) 



/(*) 



(j){x) 



5 x _i_ 1 

9 3 ' 6 



f + !(2 + x)v / T^ln 



l+Vl 11 ^ 



1-V1^5 
1 



•| - f + 1(2 + x)y/l-x arct & _ , 



for x < 1, 

for x > 1, 



-(2 + x)VT^, 
6 



X/1,T 



4 ^,r 



In the second order it is enough to take only the logarithmic term from the electron contri- 
bution 



n 2 (3) = ^(L-ivr) + C(3)-^ 



(52) 
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A discussion of the 0-meson contribution to the vacuum polarization operator, which is 



important for y/s ~ m^, can be found in |23], |24 



Appendix 2 



Here we present the direct evaluation of the collinear region contribution to the Bhabha 
scattering cross-section. Let us write the contribution of the collinear kinematics in the 
form: 



(da; 



coll 



da k\\p_ + da fcnp 4 



da fcnp'_ + da 



k\\p> 



da a + da b + da c + da d 



(53) 



For the case of photon emission along the initial electron we have 



W n = — 



uo 2 l-(3c 2 



T 



l + (l-z) 

1 — X 



-R{ Sl ,t a ,u a ) 



dr a = ^^dfi_, (54) 



uj a, 



a 



m 



c-2 



cos(fcp_), 13= \ I ;r, u = k° = xe } s l = s(l-x), t a 1 =t a (l-x) 



u 



u a (l-x), t a 



;i-x) 2 (i- c ) 



u 



tin 



-(IzfKllfO, s = iE \ 



On 



a a = 2-x(l-c), y1 = , c = cos(pTp_)- 

a„ 



Performing the angular integration over photon angles inside the narrow cone, surrounding 
the direction of the initial electron beam, we get 



Jw a dT a = 4K— dQ_ 



da; 



dci 



da; , „ yl ( _ . 9'n 



i-el/2 



i -p Cl 



4 7T ^L d Q_Hl L + ln^ )+O(6 2 ). 



x 



(55) 



We neglect the terms proportional to 0$. Collecting all the factors and reminding the con- 
tribution of the terms proportional to m 2 (see Eq. fl3~T|)), we obtain the contribution of the 
first collinear region: 



dcr„ 4a 2 a r dx 



dfi_ 



S TX J X 

A 



-4 --4 4 



X 



[ag + (l-c) 2 (l-x) 2 -q a (l-c)(l-^ 
a 2 (l-x) 2 (l-c) 2 



(56) 



For the case of photon emission inside the narrow cone, surrounding p, , in a similar way 
one gets 



dab 
dfi_ 



4a 2 a f dx 



S TX J X 
A 



-4—4^ 



X 



\a\ + (1 - c) 2 (l - xf - a b (l - c)(l - x)f 
a 2 (l-c) 2 



a b = 2-x(l + c). (57) 
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Here we used the following formulae: 

2 1 „ l + (l-x) 2 



W h 



Vi = — 



u 2 1 - (3c+ ' 
p°_ 1 — x 






1 — x 



d 3 k v b 

R{ Sl ,t b ,u b ), dr b = — ^dft_ 



lu at 



l + c)(l-x) 4b , , 



1-c 

-Si , u b = -s 1 - 



t b + M b + Si =0. 



For the cases k || p'_ and A; || p' + quantity R (if suppose 11=0) is simple 

1 /3 + c 2 \ 2 



Re. — Rd 



4 V 1-c 



dl\ rf = e 2 dx dfi_ d0i dci 

70 



xy 



2/' 



p: 



P 



2 — x + xci 
2(1 -x) 



feypi 



2 — x + xc\ 
2(1 -x 



T _ T l + (l-x) 2 
J-c — J-d — : tic 

1 — X 

C\ = cos kp' , 

1 — X, 



Cl-^-1 



^ c • (V- 



5 fellPV 

w d • ( v + 



+ Cr 



2 — x + xc\ " As 2 1 — x 
2(1 -x) 



X 



Note that kp' + = 2e 2 {\ — y). In the evaluation of the rest multipliers for these cases one has 
to be careful: 



WAT, 



l-0g/2<ci<l 



Jw d dr d 



-l+9g(l-a;) 2 /2>ci>-l 



dx _ 1 — x 

2vr— dfi_ 

x 2 



^2 

L + ln^ + 21n(l-x) 



Note that the collinear region o? is defined by the condition 1 — 6q/2 < cos kp' + < 1, which 
leads to the bounds on c\ shown above. So, the contributions of these two collinear regions 
are 



da c + da d 



a 



2 / 3 + c 2 \ 2 a }dx 



4s \ 1 - c 



IT J X 
A 



1-X + 



X 



X 



' 2 \ T 2 

L-l + ln^ + 21n(l-x) + — 



(5? 



Note that there is an asymmetry between the contributions due to the emission along the 
directions of the (initial or final) electron and the ones due to production of collinear photons 
along the positron momenta. The symmetry was broken when we decided to write a differ- 
ential cross-section with respect to the electron scattering angles (dfi_ = dcos(pTp_) d(f>). 
After an integration over a symmetrical angular acceptance the contributions would become 



equal. Compensating terms are to be extracted from the Eqs. (|56|j57|j58|) by omitting the 
terms proportional to (L — 1). 
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Appendix 3 

The following expression could be used instead of WT in Eq. (|36| ) for more precise defini- 
tion of the contribution due to non-collinear hard photon emission in large-angle Bhabha 
scattering: 

(SS) (S 1 S 1 ) (TT) 



(WT)u 



|1 - Tl(sWsx'-X' + |i - n( Sl )| 2 SlX _ x+ |1 - U(t)\H X+X ' + 



(59) 



|i-n(t 1 )| 2 * lX _x'_ 

(S3! 



+ Re 



[l -u(t))(i -u(t 1 )ytt lX -x'-x+x' + 



+ 



;i - n( s ))(i - n(s 1 ))*ss lX -x'-x + x' + 

(TS) 



(Ti5i 



l - n(t))(i - u(s)yts X '-x + x' + (i - n(t x ))(i - n( Sl ))*Mix-x'^ 

(Ti5) (T50 



where 



(1 - n(ti))(l - U(s))H 1 s X -X'-X' + (1 - n(t))(l - Yl(si))*ts lX -X+X+ J ' 

(SS) = (S 1 S 1 )=t 2 + t 2 1 +u 2 + u 2 1 , (TT) = (T 1 T 1 ) = s 2 + s 2 1 + u 2 + u 2 1 , 

(SSi) = (t 2 + tl + u 2 + u 2 1 )(t X+X ' + + t lX _ X '_-u X+X '_-u lX _ X ' + ), 
(TTi) = (s 2 + sj + u 2 + u\)(u X+X _ + u lX -x'+ + SX-X+ + S1X-X+), 

(TS) = --(u 2 + u 2 1 )[s(t + s 1 )+t(s + t 1 )-uu 1 ], 
(TS^ = --(u 2 + uj)[t(s 1 +t 1 ) + s 1 (s + t)-uu l ], 
(T X S) = -(u 2 + uj)[t 1 (s + t) + s(s 1 + t 1 )-uu 1 ], 

(TiSt) = -(u 2 + u 2 1 )[s 1 (s + t 1 )+t l (s l +t)- UUl }. 

We checked analytically that for the switched off vacuum polarization the above formula is 
equivalent to the multiplication WT in Eq. Q3TJ) : 

(^T) n |n=o = WT. (60) 

In the compensating terms, entering into Eq. ([36]) the vacuum polarization corrections, 
have to be inserted also. That can be easily done starting with Eq. (SBJ). We get 



comp 
60l 



a f dx 



TT J x 

A v x 

dgb(MU 

x dfi_ X 



1 — x H In 

2 i 



<9g(l -x) 2 x 2 



1-x) 



+ 



'l-x + ^ln^ + ? 



2i 



X 



dcr (l-X, 1) dcr (l,l -x) 

dfi_ 
dor (l,l)8a 



dfi_ 



, , 6» Ae 1# _ 
,o ln ( ct go) m — dfi_. 

dii_ 7T 2 £ 



(61) 
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We note that quantity R (see Eq.(^)) in the ultra-relativistic limit s ^> m^ can be derived 
from Eg. ((59]) . One has to omit there all terms except the ones proportional to (SS), (SS\), 
(Si Si) and divide by 4. 
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Figure 1: The differential forward-backward asymmetry for e + e — > yU + yU . The parameters 
are e = 0.51 GeV, ^ = 10°, e th = 0.2s. 



Figure 2: The differential forward-backward asymmetry for e + e — > yU + /i . Beam energy 
e = 1.55 GeV, other parameters the same as in Fig. 1. 
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